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We consider the evolution of cosmological perturbations during radion stabilization, which we 
assume to happen after a period of inflation in the early universe. Concentrating on the Randall- 
Sundrum brane world scenario, we find that if matter is present both on the positive and negative 
tension branes, the coupling of the radion to matter fields could have significant impact on the evo- 
lution of the curvature perturbation and on the production of entropy perturbations. We investigate 
both the case of a long-lived and short-lived radion and outline similarities and differences to the 
curvaton scenario. 
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I. INTRODUCTION 

Present and future cosmological observations enable us to constrain or rule out models of the early universe. For 
example, the observations made by the Wilkinson Microwave Anisotropy Probe (WMAP) already put some constraints 
on inflationary cosmology, 0, S S 3 ■ The observations are consistent with purely adiabatic perturbations with a 
scale invariant power spectrum, but some "anomalies" , such as the small power of the quadrapole and features in the 
CMB anisotropy spectrum, have been reported, 1^. 

Future cosmological observations however, will not only constrain current popular models of the early universe, they 
will also test models "beyond the standard lore" , such as models based on string theory, extra dimensions, etc. For 
example, alternative mechanisms to inflation have been proposed recently. In contrast to the standard inflationary 
picture, where density perturbations are generated during inflation and then stretched onto superhorizon scales, flelds 
other than the inflaton field generate the initial perturbations in the new mechanism. In the curvaton scenario initial 
isocurvature perturbations are transformed into curvature perturbations by a subsequent decay of a second field. This 
field, dubbed the curvaton, is already present in the early universe but is dynamically uniinportant during inflation, 
■ leHvHiSmillllllll. Another example is the idea of "modulated perturbations" , E 111 [jj, IjJ, [jj . According 
I to this idea, coupling constants (and other physical properties such as masses of the particles) are functions of the 
vacuum expectation value of fundamental (light) scalar fields present during infiation. Because of vacuum fiuctuations 
in these fields, different regions in spacetime have different values of coupling constants. These fluctuations can be 
converted into curvature perturbations, for example during reheating which is not homogeneous in space if the decay 
rate itself fluctuates. Note that both the curvaton scenario as well as the idea of modulated perturbations need 
additional scalar fields in the theory, which are dynamically unimportant- at least initially, during inflation. 

In this paper we do not seek an alternative to the standard picture of how perturbations are generated, but rather 
investigate how far the stabilization of moduli fields can alter the evolution of cosmological perturbations generated 
during a period of infiation in the early universe. Moduli fields appear in theories beyond the standard model based 
on supersymmetry or superstring theory. In some models they can be long-lived and the life-time can be even larger 
than the age of the universe. The prime example of a moduli field is the radion, measuring the size of the extra 
dimensional spacetime in brane world scenarios (see e.g. pol l2ll for recent reviews on brane cosmology). One 
important property of the radion is that it couples explicitly to matter fields on the branes (see e.g. |23j and references 
therein). Depending on the details of the bulk geometry and the matter content on the branes, the coupling is field 
dependent and can become quite large. It is the coupling between the different matter forms to the radion which 
can lead to interesting consequences for the evolution of perturbations. Furthermore, the radion can be short or 
long-lived depending on the stabilization mechanism. In the case of the Goldberger-Wise stabilization mechanism in 
the Randall-Sundrum brane world, ,2^j, the radion obtains a mass of the order of TeV, which makes is short-lived, 
|25|. In this paper we investigate the impact of radion stabilization on the cosmology and dynamics of cosmological 
perturbations. We consider both the case of a short-lived and long-lived radion; in the latter, the radion will constitute 
some (or all) of the dark matter in the universe. Throughout the paper we work in the Randall-Sundrum scenario, 
|26j. because the coupling of the radion to matter on both branes can be quite large and new effects will appear (see 
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e.g. [2?, for the effects of the radion on the CMB anisotropies) . In p^lsoHsif the radion and its role in cosmology 
were discussed in considerable depth. Here we take into account the warped geometry of the bulk. 

The paper is organized as follows: in Section we will present and describe the low-energy effective theory of the 
Randall-Simdrum brane world. We then present the equations of motion for both background and perturbations. 
In Section IIIII we will discuss the evolution of the radion and the different constraints on the theory coming from 
nucleosynthesis, overclosure of the universe and the amount of entropy perturbations generated. In Section llVI we 
present our numerical results of the evolution of the curvature and entropy perturbations. Our conclusions are 
presented in Section Ivl 



II. THE EFFECTIVE ACTION AND EQUATIONS OF MOTION 

The form of the low-energy action for the Randall- Sundrum brane world scenario, in which two boundary branes 
are embedded in a slice of an Anti-de Sitter (AdS) spacetime, has been studied in the literature already, see e.g. 
[2^ 112, 133L I33 . '3^ and references therein. At low energies the theory is a scalar-tensor theory with specific matter 
couplings. In the Einstein frame it takes the form 



^EF = d, xy/—g 



1 .n~hdRf-v{R) 



IGttG 2 

+Si^H^i,A'{R)g,,) + (*2, B'{R)g,,). (1) 

We have allowed for the possibility of two matter forms, each of which is confined to the positive or negative tension 
brane. The couplings between the fields on the branes and the field R are described by the functions A{R) and 
B{R). The coupling to matter has some interesting consequences, the primary feature being that we no longer have 
energy-momentum conservation for each component. In fact 

"^t^Tt" = a'R{d''R)T,, (2) 

where i describes the particular brane in hand. Note that the field couples to the trace of the energy-momentum 
tensor and so there is no coupling to radiation. One defines the coupling functions 

(1) d\nA (2) d\nB 
= = -dR- 

The energy conservation equation (0) is modified in this theory because the masses of particles on the branes are 
functions of the radion, and hence can vary with time if the radion varies with time. Alternatively, we could choose 
a frame in which the masses of the particles on the positive tension brane are constant. In this frame, the masses of 
the particles on the negative tension brane as well as the four-dimensional Planck mass vary with time. We will study 
the theory in the Einstein frame, because after stabilization of the radion the frames agree. 

One problem common to nearly all models with moduli fields is how to stabilize these fields. We shall include a 
potential of the form 

V{R) = ]^Ml{R-R,)\ (4) 

but here we note that the origin of such a potential might have to be derived from non-perturbative effects in the 
underlying theory. As we will see below, the mass Mji is constrained by requiring that the field does not overdose 
the universe and also by nucleosynthesis. 



Finally, the couplings A{R) and B{R) have the form 





"fl^ = ^ coth I I . (6) 



The geometry of the higher-dimensional spacetime is encoded in these coupling functions. Depending on the value 
of i?, the couplings can be quite large. For small values of i?, the function a\^' behaves as Up,' ~ whereas one 
finds a^^^ - R/M^i- 
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The theory described above is a good description of the two-branc system as long as the energy densities on both 
branes are much less than the brane tension. All length scales considered should be much larger than the curvature 
scale of the Anti-de Sitter bulk. Furthermore, we have assumed that Kaluza-Klein modes are irrelevant at the energy 
scales we consider. 

Primarily we are concerned with the application to cosmology. We allow matter to be present on both branes, each 
with an energy density pm'' and p"m respectively. We also include radiation with energy density p^. To end up with a 
realistic cosmology, i.e. with late time acceleration, we include a cosmological constant A. Assuming a flat universe, 
the equations for the background read: 



1 



PA 



i? + (3i/ + T)R 



^+C^(l)«(l)+a(2).. 



pW = ^3Hpl:^ +o^Rp(^, i = l,2, 



-AHp^+TR^. 



(7) 

(8) 

(9) 
(10) 



We will define the density parameter for the different components by 17^ = Pi/ Per as is usual. 

The equations of motion of perturbations in matter, the radion field and the metric can be obtained from the 
effective action. We consider scalar perturbations only in this paper and work in the longitudinal gauge, in which the 
perturbed metric has the form (neglecting any anisotropic stress). 



ds'^ = -(1 + 2^)dt^ + a'^{t){l - 2*)rfx2. 
The energy-momentum tensor for the two fluids has the form 

Tf^ _ I -{P + ^P) --^d'Sq 



(11) 



(12) 



According to our assumptions, Spi = LOiSpi for both matter and radiation. Then, expanding in Fourier modes, the 
perturbation equations read at first order 



6R + {3H + T)dR + -^6R ■ 



SR 



-a 



R "Pm 



a 



R "Pm 



(1) 

dR ^"^ 




Sp<;i^+SH6p: 



U) _ 



3*p« + 



-2* — 

\dR 



i = l,2 



Spy + AHSp-y - 4*0^ + -^Sq^ 



Sqii^+3HSqii^ 
5q-y + SHSq^ 



'^p 



+ -^^Pi 



3if(* + il*) - 

a"' 



* = - 



Qj^ + SQ-y, 

-urp^^ 5R, 
0, 

1 



Sp. 



where we have used the abbreviations 



a'i'5R + 



da 



OR 



■RSR 



6Qy 
Sp 



TR^, 

2TR5R - *r.R2 

bpy + dp^^^ 



pfi + o^m(Ji 



dV , 



0) A 

m ' J 



1,2. 



dp'i^ + RSR + ^SR - ^R'^. 



(13) 
(14) 

(15) 
(16) 

(17) 

(18) 

(19) 

(20) 

(21) 
(22) 

(23) 
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The fluids and the radion field exchange energy during the cosmological expansion due to the non- vanishing coupling. 
This is to be compared with the original curvaton scenario, in which the curvaton field decays into radiation (and 
matter) . 

Furthermore, there is also energy transfer between the radion field and the radiation fluid due to the F interaction 
term. The decay rate is given as, |36|. 

F = Q , , (24) 

where the constant a is determined by the number of extra dimensions and the compactification and is 0(1). We 
shall set a = 1 for the duration of this paper. This means the radion has a decay time of 



(25) 



Thus for Mr < 0.1 GeV, the radion would effectively be stable over the lifetime of the universe. If we require the 
radion to decay before nucleosynthesis, then we need 

r>i/(t„uci)-10-^=^Mp,. (26) 

This would then require a radion mass ^ 5 x 10"^** Mp; « 10^ GeV, assuming the constant a is of order unity. 
The radion stabilizes at ~ so this gives a long period of radion domination. 

The evolution of the perturbation quantities is described in full by equations H13I18|I . However, it will be helpful 
for our understanding to recast these. In general, the density perturbations, 5pi, and the metric perturbation, 4*, are 
gauge dependent. It is more common to work with gauge-invariant quantities. One important quantity, the curvature 
perturbation on constant density hypersurfaces, is (see e.g. |37L ISSL l39| ) 

^ = 5r + i/^. (27) 
P 

Note that this differs in sign from some of the definitions in the literature. In addition we can define the curvature 
perturbation on uniform j-fluid density hypersurfaces as, [s^ . 

C, = 1' + i/^. (28) 

P^ 

With this prescription, it is clear that the total curvature perturbation is a weighted sum of the individual components, 

C = E-C- (29) 

i ^ 

Furthermore, we are able to define a relative entropy perturbation between two components as 

5., =3(C.-0)- (30) 

To make some predictions about the behavior of the perturbations we work in the "separate universes" picture, 
which effectively allows us to ignore the spatial gradient terms. In the long-wavelength limit, the perturbed continuity 
equation becomes 

5p + iH{5p + 5P)^i{p + P)^, (31) 

where we define 

P = Ep„ P = Ep», <5p = E(5p„ 5P = Y,5P^. (32) 

i i i i 

If one rewrites equation H31(l in terms of the total curvature perturbation, ^, it is relatively straightforward to show 

C= ^IJ-^Pnad, (33) 
P + P 
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where the non-adiabatic pressure perturbation is SPnad = SP — c^Sp and the adiabatic sound speed is = P/p- This 
means that the total curvature perturbation will be constant on large scales for purely adiabatic perturbations. In 
the context of brane worlds, where energy can be exchanged with the bulk space time, this quantity is not necessarily 
constant (see e.g. "i^, "4^, "i^ ) . 

It is possible to break this down further so that one may write, p^. 

JFnad = SPnt+5P,ch (34) 

= ^(<5F,-c2<5p,) + ^E'^»'^'j(c?-s')'5^.' (35) 

i i,j 

where cf = Pi/pi is the adiabatic sound speed of the fluid and is related to Cg by 

= E J^- (36) 

i 

Further, let us define 

5PinM = - C^Jp,, (37) 

SO that SPint = J2t ^Pint,i- 

On superhorizon scales it is possible to write the individual perturbation equations in the form, [iflj. 



Ci = — i^Pi - Ci 6 Pi) H — 

Pi Pi 



SQ^ I P Qi \ Sp, ^ 

Qi [2p Q^) Pr H 



The Qi and 5Qi are given in equations (|19I22|) with Qr = —Q')n — Q^^m — Q-y and similarly for SQr. From this, one 
can immediately see that for non-interacting perfect fluids with Qi = and 5Pi = cfSpi, the individual curvature 
perturbations for each fluid remain constant on large scales. 

There is now large scope for the setting of the initial conditions of the perturbations. We shall assume that the 
initial conditions result from a period of inflation in the early universe. Then, it is known that any light scalar field, 
will pick up a perturbation, [40ll4l|. 



(38) 



0, (39) 



where this is evaluated at horizon crossing for the inflaton field. This depends on the scale of inflation but if this is 
to be taken as V < 1.5 x 10"^ M^j, 0, this tells us that 

|(5i?| < 1.1 X (40) 

If this results in too much entropy production, this will offer a constraint on the scale of inflation. The remaining 
initial conditions are arranged to be adiabatic initially, as given in 45] , 

* = Spii^ = -2Cpl^J, 5p^^-^Cp^, (41) 

where C is some initial scale to be decided in advance. We shall set C = 10^^ in accordance with current observations. 
Since we are dealing with linear perturbation theory, it is only the ratio \6R\/C' that is important rather than their 
absolute values. In fact, in future we shall normalize the curvature perturbation and plot Ci/Cini- Note that these 
initial conditions are different from those chosen in where 6R = initially. For the duration of this paper we shall 
set 

(5i? = -1.1 X 10"^ Mp/, (42) 

since this will give the largest possible effects. 

With this prescription, it is worth discussing what we might expect to see. From equations and one 
can see that ^ is sourced in two ways. Firstly when Pi ^ Pi(pi), which makes the flrst term in H35|) non-zero, and 
secondly if there is an entropy perturbation present. For perfect fluids, (5Pintr,i will be zero but for a scalar field this 
is not necessarily the case. Therefore we will always generate evolution in (^n and thus in ^. Furthermore, our initial 
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conditions, whilst adiabatic between the fluid components, necessarily produce an entropy perturbation between the 
radion and the fluids. This all holds without any couplings. One should add that since the total curvature perturbation 
is a weighted sum of the individual components, the efl'ect of the scalar field on ( is very much dependent on its energy 
density fraction. During radiation domination, which will be constant- if F 7^ this no longer holds. During 

matter domination it is then described by 

C«f^We+f^(^)e)+17HCH. (43) 

As one moves into the A dominated phase this should then become constant again. If couplings are introduced then, 
from equation (|38|l . we know that all the d evolve. During the energy transfer between the fields, (5Pint,i 7^ and we 
will also produce entropy perturbations. We should expect to see further evolution in which will be dependent on 
the size of the coupling terms. 

If one includes coupling between the radion and radiation only (F ^ 0,a^^ = 0), we should see an effect in (^r, 
and therefore C, as soon as the radion stabilizes, because it becomes the dominant component. All of its energy is 
transfered to radiation as it decays. Once it has decayed, gives no contribution to the total. However, since there 

(i) . 

is no coupling to matter, the have remained constant and C will return to this value. Of course, this conclusion 
will be modified if the coupling terms are non-zero. 

III. COSMOLOGICAL EVOLUTION OF THE RADION 

Before we discuss the evolution of the perturbations, we will consider the evolution of the radion during the infla- 
tionary epoch and the subsequent radiation and matter dominated epochs. We also discuss cosmological constraints 
the theory has to fulflU. We begin with the inflationary era. 

A. The Epoch of Inflation 

For definitiveness we consider an inflaton fleld conflned to the positive tension brane and a chaotic inflation potential 
of the form V{(f>) = M|0^/2, where is the mass of the inflaton fleld. Further, because we want to study purely 
adiabatic initial conditions between matter and radiation, we assume that during the last 60 e-folds the inflaton fleld 
dominated the energy density of the universe. We are not interested in the case in which the radion was the inflaton 
fleld. With these assumptions, there are three cases of interest. In the flrst case, the mass of the radion is much larger 
than the mass of the inflaton fleld, i.e. Mr ^ M^. Depending on the initial conditions for R, the radion fleld could 
have dominated the initial expansion of the early universe. In the second case, the masses of the radion and inflaton 
fleld are comparable. In this instance, the energy density of the radion must have been much smaller than the energy 
density of the inflaton to comply with our assumptions. Finally, the radion can be much lighter than the inflaton 
fleld. Let us briefly discuss the three cases separately. 

1. Mr < 

In this regime, the radion is very light compared to the inflaton fleld. Its potential energy can be neglected and the 
radion fleld is driven by the coupling to the inflaton fleld only. As established earlier, |46l |. the radion fleld is driven 
towards small values. 

2. Mr « M^ 

Here a period of double inflation can be realized, where each phase is driven by either the radion or the inflaton 
fleld. In accordance with our assumptions, we consider the case in which the second period of inflation, driven by the 
inflaton fleld, is longer than 60 e-folds. With this prescription, the radion and the inflaton fleld stabilize at the same 
time. 

3. Mr > M^ 

This case is very similar to the second example, except that the radion now stabilizes much earlier. Initially, it is 
possible that the radion dominated the expansion rate of the universe before settling into the minimum of the effective 



7 



potential. In principle, the radion could have set up the initial conditions for a period of inflation. 

In the following, we consider a scenario where the radion decays well after inflation, which corresponds to the first 
case only. This ensures that the perturbations created during inflation are adiabatic as is the case in single field 
inflation. 



B. The radiation and matter dominated epochs 



The behavior of the radion field depends on its mass. If it is very light, the lifetime can be very large and in 
fact, in this instance, we assume that the lifetime is much larger than the age of the universe. However, depending 
on the origin of the potential for the radion, it could also be very heavy. In the case of radion stabilization via the 
Goldberger-Wise mechanism the radion mass is of order TeV. Here, the radion decays into radiation and matter in the 
very early universe. We will consider these two cases in the rest of the paper. Further, we need to make assumptions 
about where the matter fields are located. In principle, both the standard model particles as well as the dark matter 
particles can be located only on the positive tension brane, whereas there is no matter on the negative tension brane. 
Alternatively, on can imagine the situation where only the standard model particles are confined on the positive 
tension brane, whereas the dark matter particles are confined on the negative tension brane. Although mixtures of 
these scenarios can be imagined, in this paper we will allow only baryons to live on the positive tension brane; the 
dark matter being made up from matter on the negative tension brane and from the oscillating radion field. 

In order to make our study compatible with present day observation, we take values for the parameters as given by 
the most recent WMAP data, |^. This means that today our universe should be almost flat and made up of 

r^A = 0.73, = 0.23, = 0.04. (44) 

Furthermore, we ensure that matter-radiation equality occurs in the range obtained by the WMAP team, 

3200 < zeq < 3400. (45) 

Nucleosynthesis makes excellent predictions for the abundances of the light elements and we require that these 
predictions are not modified too much. Nucleosynthesis gives a constraint on the expansion rate -ffnuci at the time 
when the light elements form. If the particle masses change during the history of the universe, as is the case in our 
theory, this will infiuence the time dependence of H and therefore the expansion rate at the time of nucleosynthesis. 
The evolution of the masses is specified by the functions A and B. The requirement that any changes of H^^ci are 
constrained by nucleosynthesis gives, [i^ . 

Anucl - ^0 ^ ^^g^ 



^0 

where Anuci and Aq are the coupling functions at the time of nucleosynthesis and today respectively. A similar bound 
is valid for B if some matter form is present on the negative tension brane. Note that these considerations neglect 
variations in the reaction rates, which itself are functions of masses- and therefore of A- if the standard model particles 
live on the positive tension brane. In contrast, if they live on the negative tension brane instead, the reaction rates 
would be functions of B. Therefore, the constraint above is a rather conservative one. 

If the radion is light [Mfi < 0.1 GeV), we have seen that it is effectively stable over the lifetime of the universe. It 
can be shown, j42|], that for a scalar field, ip, oscillating in the minimum of its potential, V ^ ip^, the energy density 
scales 

PV---a"^"/'"+'^. (47) 

In our case this would mean the radion, R, behaving like non-relativistic pressureless matter. However, we also have 
the extra source terms, a^^^ p^}} , on the right hand side of equation ((HJ that could amend this result. The oscillations 

of R are damped by the Hubble parameter and the energy density p^m will scale approximately like . This means 
that the source terms will decay faster than the others in equation ||SJ). Therefore, we should still find that the 
radion behaves like matter when it oscillates in the minimum of its potential, regardless of the new couplings. This 
is confirmed by our numerical results. In fact, the couplings are most important for the background solutions as the 
radion rolls down the potential. In practice, this transfers energy between the different species. 

Since R behaves like matter when oscillating at its minimum, this allows us to investigate a number of possibilities: 



1. R plays the role of dark matter alone. 
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2. R is part of the dark matter with the remaining contribution from the positive and/or negative tension brane. 

3. R is subdominant with the matter on the second brane making up aU of the dark matter. 

Stabilization occurs when H ^ Mji and from this point pn ^ (which inevitably succeeds ~ a""*). Therefore, 
one can not stabilize the radion too early as this would lead to premature matter domination. This means that the 
radion cannot be too heavy. Moreover, if the radion is very light, then the background dynamics are more sensitive 
to the couplings. This may result in too large an evolution of R. Some of this can be compensated for in the initial 
conditions, but in general we shall find it difficult to generate extremes of masses. If R is to play the role of dark 
matter, then its mass should be around Mr ~ 10^^"Afp/, ^^|Eq|> ^o agree with galactic halo results. We shall take 
this to be given in this instance and match our initial conditions to the correct cosmology. The situation in which R 
is always subdominant may be interesting because even though we can not see the field, its coupling to matter may 
have some interesting effects on the perturbations. 

Whilst it is possible to have dark matter- in addition to baryons- present on the positive tension brane, the coupling 
function in equation Q is small in general. In contrast, it is possible to produce large couplings on the negative tension 
brane, as it can be seen from equation ©. 

IV. RADION STABILIZATION 

Having studied the background evolution, we now turn our attention to the evolution of the perturbations. We first 
discuss the case of a long-lived, i.e. stable radion field, before we consider the effects of radion decay into radiation. 

1. The case of a long-lived radion 

In the Randall- Sundrum model, the coupling of the radion to matter on the positive tension brane is never large, 
c^n^ < l/\/6- In contrast, the coupling to matter on the negative tension brane has the property a^'' ~ 1/R for small 
R. This allows the possibility of generating rather large couplings. In our first example, we take the radion to be 
approximately half of the dark matter. The results are shown in Figure ^ Here we are able to generate couplings 

(2) 

as large as a)^' ^ 1000 and unsurprisingly, this gives some strong effects on ^. One sees that there are large entropy 
perturbations and some interesting changes in the individual Ci ■ The most striking feature is that C evolves to negative 
values. 

Let us now consider a subdominant radion and examine the entropy production from the couplings. The results 
for this setup are shown in Figure 13 where we stabilize the radion at i?c = 1 Mpi . The couplings can have important 
consequences for the background evolution if the radion is very light. These couplings dictate the evolution of the 
radion initially and push the radion up its potential. This can lead to a period of radion domination as it gains energy 
from matter. We find from our numerical results that if R is to remain subdominant we require AIr > 10^^^ M^;. 
Once more we generate entropy perturbations. 

We have seen that for a long-lived radion its stabilization has some interesting consequences for the history of the 
universe. Even if it remains subdominant in the background evolution, it is still possible for it to generate entropy 
perturbations between matter and radiation due to the couplings. 

2. The case of a short-lived radion 

We move on to discuss the case of the decaying radion. In order for this to occur before nucleosynthesis- so that the 
universe is radiation dominated- we require Mr > 10-"Mp,. This should give us a large period of radion domination 
which is followed by radion decay into radiation. In the case of the Goldberger-Wise stabilization mechanism via a 
bulk scalar field, the mass of the radion is of order 1 TeV and therefore decays naturally in the early universe. In 
some sense, it could be a candidate for the curvaton. Here, however, we are concerned how initial perturbations from 
inflation are modified by the stabilization and decay of the radion. We will set the mass of order 100 TeV so that the 
decay occurs before nucleosynthesis. 

Again, the background initial conditions are set to produce the universe that we observe today. As we have already 
discussed, C should evolve when R becomes dominant and then again at matter-radiation equality. If a]^ = 0, then ^ 
should return to unity. Any deviation from this is down to the couplings to matter. This can be seen from equation 
(|43|l because flR — and (^m are constant. 
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FIG. 1: This plot shows the results for a universe with Qa = 0.73, Qr = 0.13, = 0.10, Sl^' = 0.04. Furthermore, we 
set Mr = y/6 X 10"^*^ Mpi . In the top left we show the evolution of the density parameters, Qi. In the top right plot we see 
the evolution of the total curvature perturbation, which demonstrates the effect of the couplings Q^', compared with the 
uncoupled case. We set Ro = 0.10 Mpi and = 0.03 Mpi so that 10"* < a^^' Mpi < 10~^ and 10 < a^^' Mpi < 10^. The lower 
left plot shows the behavior of all the individual (^i and on the right we show the evolution of all entropy perturbations, Sij, 
between the different species. 



In Figure 131 we stabilize the radion close to the horizon of the AdS spacetime at i? = to generate large couplings 
on the negative tension brane. One can see that the background dynamics are as expected and that the radion decays 
quickly into matter before nucleosynthesis- this is ensured by setting Mr — V6 x lO^^^Mp;. The are sourced 
when the radion is stabilized and then remain constant, whereas is sourced during the decay of R. This is reflected 
in the not insignificant production of entropy. 

By means of comparison, we show the equivalent plots in Figure 0] for no coupling to matter. Whilst the picture is 
similar, we find ( returning to unity at matter-radiation equality. Furthermore, there is less entropy production than 
before. 



V. SUMMARY AND DISCUSSION 



In this paper we have studied the evolution of cosmological perturbations during radion stabilization. We have seen 
that stabilizing the radion field, which describes the inter-brane distance, can have non-trivial effects on the curvature 
and entropy perturbations. We have studied two situations, a long and short-lived radion. In the short-lived case, 
we allow the radion to decay into radiation only with the radion coupling to matter naturally. In setting the initial 
conditions for the perturbations, we assume that our universe has already undergone a period of inflation. The scale 
of inflation determines the initial conditions for the radion- in particular the amplitude of the field fluctuation- with 
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FIG. 2: In this example, we generate a universe with SIa = 0.73, Qr — 0.0001, fim' = 0.23, flm^ — 0.04 today. We have set 
Mr = \/6 X lO"^'' Mpi. The difference between this and the previous FigureQis that the radion field is subdominant throughout 
the history of the universe. This is clear from the top left plot in which we show the evolution of the density parameters, Qi. 
In the top right plot we show the evolution of the total curvature perturbation, ^, with and without coupling terms. We set 
Ro = 1.002 Mpi and Rc ~ l.OAIpi so that a^-* Mpi ~ 0.15 and a^' Mpi ~ 1.05. As one can see, the evolution of the curvature 
perturbation is significantly modified, even when the radion field is subdominant. The lower left plot shows the behavior of the 
individual (i and on the right we show the evolution of the entropy perturbations, Sij. 



the fluid perturbations set to be adiabatic initially. This will naturally source the total curvature perturbation, since 
there is an initial entropy between the fields and the fiuids. Additionally, energy is transfered between the radion field 
and matter on each brane due to the coupling terms. This enhances the amount of entropy produced. The amount 
of entropy production coming from the coupling terms depends on the initial conditions of the radion field just after 
inflation as well as on the field value at stabilization. 

The evolution of ^ and the Q- and therefore Sij- are influenced by a number of parameters. The values Rq and Rc 

are important as they determine the size of the coupling terms cx^j^ . In addition, |(5i?ini| and C determine the initial 
entropy perturbation between radiation and the field which is transfered into part of the curvature perturbation. This 
is similar to the curvaton scenario. The difference here is that there is an initial curvature perturbation and even if 
the field decays into radiation only, there is still an energy transfer to matter via the couplings. 

The constraints on the theory come mainly from the CMB anisotropics. The effect of the radion on the CMB 
anisotropics has been investigated earlier, |27l l28|. but in this paper we have included an initial field perturbation, 
SR ^ 0, set by the infiationary scale. The large-scale temperature fluctuation can be written as, [sif. 

+ + I^S^,] , (48) 

PM ' PM ' PM I 




since the radion behaves like matter at matter-radiation equality in all cases we have considered. In this expression we 
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FIG. 3: In this instance, we consider a short-lived radion field with a mass Mr^VGx lO'^^Mpi. We set Ro = O.lOA'/pi and 
Rc = 0.03 Mpi so that 10"^ < a^^' Mpi < 10"^ and 10 < a^^' Mpi < 10^ This results in a universe with = 0.73, = 0.23, 
^rn ~ 0.04. In the top left we show the time dependence of the 57^. We see that after a period of radion domination it decays 
into radiation. In the top right plot we see that evolution of the total curvature perturbation, The lower left panel shows 
the behavior of the individual ("i and on the right we show the evolution of all of the entropy perturbations, Sij, between the 
different species. 



have set pM = Pm + Pm + Pr and we neglect any neutrino contribution. Observationally constrained is the "effective 
baryon isocurvature perturbation" , defined as 

(2)c(2) I c 

OcB = J;,i4 H TjN , (49) 

Pm 

where in this expression we have assumed that the matter on the first branc is comprised of baryons only. Since the 

(2) 

amplitude of Smy as well as of Sr^ depend on the details of the stabilization mechanism, such as sizes of the couplings 
during radion evolution and the mass of the radion, the CMB constrains these parameters to a certain extent, in a 
similar manner to the curvaton scenario, .51;]. 

It is possible to extend our work in several different ways. Firstly, the initial conditions between the matter and 
radiation do not have to be necessarily adiabatic. In fact, the radion is a candidate for a curvaton field, if inflation 
happens at a higher energy scale than the energy scale of the radion potential. In this case, the radion field picks 
up perturbations during inflation which can be transfered to a curvature perturbation during radion stabilization. 
Since in the Goldberger-Wise mechanism the mass of the radion is of order TeV, this is a realistic scenario to be 
investigated. A detailed analysis would give constraints on the mass of the radion as well as on initial conditions of 
the field just after inflation. We note, that our setup is similar to the one discussed in Q. However, in our case the 

(i) 

radion field decays into radiation only and is coupled to matter via the functions a)j only. 
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FIG. 4: We make the equivalent run to Figure El but without couphng to matter. In the upper plot, we see the evolution of 
the total curvature perturbation, ^, with no couplings and the equivalent from Figure |2l The effects of the couplings on ^ play 
a role only well after radion decay. However, if one compares the two lower left plots, it can be seen that the individual C^m' 
evolve at radion stabilization. With no couplings, it is only C^-y that undergoes any change. In the lower right panel, we show 
the evolution of the entropy perturbations. 

A second point to be investigated is the effect of the potential on the stabihzation and on the evolution of per- 
turbations. If the radion undergoes minimal evolution then any potential with a local minimum can usually be 
approximated by the expansion 



V{R)^Vo + -MUR 



Rc 



(50) 



In this instance our results will should not change substantially. However, if the radion evolves sufficiently such that 
the potential can no longer be approximated by the expansion in equation (|5()|l . then the energy transfer between the 
radion and matter can be significantly different. 

Our work also has implications for other models such as those with constant couplings. If these couplings are 
large enough, the energy transfer between the modulus and matter can be large and an entropy perturbation will be 
produced, even if the field is subdominant throughout the history of the universe. 
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